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ABSTRACT 

We study numerically the anisotropy of the cosmic ray (CR) flux emitted by a single source calculat¬ 
ing the trajectories of individual CRs. We show that the contribution of a single source to the observed 
anisotropy is determined solely by the fraction the source contributes to the total CR intensity, its age 
and its distance, and does not depend on the CR energy at late times. Therefore the observation of a 
constant dipole anisotropy indicates that a single source dominates the CR flux in the corresponding 
energy range. A natural explanation for the plateau between 2-20 TeV observed in the CR anisotropy 
is thus the presence of a single, nearby source. For the source age of 2Myr, as suggested by the 
explanation of the antiproton and positron data from PAMELA and AMS-02 through a local source, 
we determine the source distance as ^ 200 pc. Combined with the contribution of the global CR sea 
calculated in the escape model, we can explain qualitatively the data for the dipole anisotropy. Our 
results suggest that the assumption of a smooth CR source distribution should be abandoned between 
~ 200 GeV and 1 PeV. 

Subject headings: High energy cosmic rays, cosmic ray anisotropies, galactic magnetic field. 


1. INTRODUCTION 

The observed intensity of cosmic rays (CR) is charac¬ 
terised by a large degree of isotropy up to the highest 
energies. The measured dipole anisotropy S increases up 
E ^ 2 TeV, forms an approximately energy independent 
plateau between E ^ 2-20 TeV, before the amplitude 
decreases again. Finally, at energies E > 100 TeV, the 
anisotropy grows fast with energy. Such a behavior of 
the dipole anisotropy 5 as function of energy seems at 
first sight difficult to reconcile with diffusive CR prop¬ 
agation; In this picture, the scattering of CRs on inho¬ 
mogeneities of the Galactic magnetic field (GMF) con¬ 
verts their trajectories to a random walk, erasing thereby 
most directional information. Cosmic rays are most ef¬ 
fectively scattered by those turbulent field modes whose 
wave-length equal their Larmor radii. Therefore the scat¬ 
tering rate is energy dependent and determined by the 
fluctuation spectrum of the turbulent GMF. As result, 
both the diffusion coefficient and the GR anisotropy are 
expected to increase with energy. This picture is sup¬ 
ported e.g. by recent data from the AMS-02 experiment 
for the B/G ratio which are consistent with the simple 
power-law D{E) oc E^-^^ up to the rigidity 1.8 TV. 

In the diffusion approximation, Tick’s law is valid and 
the net CR current j{E) is determined by the gradient of 
the CR number density n{E) = N/(Ey) and the diffusion 
tensor Dab{E) as ja = —Dab^tn. The dipole vector 5 of 
the CR intensity / = c/(47r)n follows then as 

r ^ ja \ 


Within this approximation, the diffusion tensor Dab is an 
external input. The traditional approach to determine 
Dab as a function of the assumed magnetic field uses 
kinetic theory, see e.g. [Berezinskii et akl (119901) . Such an 
analytical approach allows one to connect both Dab and 
Sa to the spectrum of magnetic field fluctuations, involves 


however the use of approximations. 

An alternative method uses the trajectories of indi¬ 
vidual CRs calculated numerically solving the Lorentz 
equation in turbulent and regular magnetic fields. This 
approach is numerically expensive and has been re¬ 
stricted mainly to the calculation o f the diffusion ten¬ 
sor Dab (iGiacalone fc Jokipiil Il999t iCasse etld] 120021 : 
iGiacinti et al.l l2012af) . Since the anisotropy in the as- 
trophysically interesting cases is small, its calculation re¬ 
quires a much larger number of CR trajectories. There¬ 
fore the determination of the CR anisotropy expected 
from a smooth source distribution has not been feasible 
using this method directly. Instea d, the approximation 
proposed bv iKarakula et al.1 (119721) has been widely used 
in the high-energy range, E > 10^® eV. In this letter, we 
concentrate on the case of a single source and show that 
the calculation of the anisotropy from first principles is 
possible for an interesting range of parameters. 

The propagation of CRs can be divided into three dif¬ 
ferent regimes: The free streaming, the anisotropic and 
the quasi-Gaussian diffusion of CRs. In the first regime, 
the change of the particle momentum is small. The inter- 
m ediate regime o f anisotropic d iffusion was explored first 
bv iGiacinti et al.l (|2012aL[201,‘lD : a study of the anisotropy 
in this regime will be presented elsewhere. Here, we con¬ 
centrate on the case of quasi-Gaussian diffusion and show 
that then the dipole anisotropy of a single source with 
age T and distance R is given by d = 3i?/(2cT), inde¬ 
pendent of the regular and turbulent magnetic field. As 
an application of this result, we consider the observed 
CR dipole anisotropy. We propose as explanation for 
the plateau in the observed anisotropy between 2-20 TeV 
that a single nearby source dominates the CR flux in 
this energy range. At lower energies, the anisotropy de¬ 
creases because the fraction of CRs contributed by the 
single source goes down. We show that the anisotropy 
data between 200 GeV and 4 PeV can be naturally ex¬ 
plained by a local GR source with the age of 2 Myr and 
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distance ^ 200 pc, consistent wi th the characteri s tics o f 
the single source determined in iKachelriefi et alJ l)2015f) 
for an explanation of the antiproton and positron data 
from PAMELA and AMS-02. 

2. ANISOTROPY IN THE QUASI-GAUSSIAN 
REGIME 

We avoid the limitations of the diffusion approxima¬ 
tion by calculating the trajectories of individual CRs in 
a given regular and turbule nt magnetic field usiiig the nu - 
merical code described bv iGiacinti et alJ (j201lL I^12bl l. 
We use as spectrum V{k) of the magnetic field fluctua¬ 
tions isotropic Kolmogorov turbulence, V{k) oc with 
a = 5/3 and set the maximum length equal to 25 pc. 
For the numerical calculations, we use nested grids. This 
allows us to choose an effective I'^i^ sufficiently small 
compared to the Larmor radius considered. 

We model the GR source as an instantaneous injec¬ 
tion of CRs at a single point. Since we are interested in 
the generic behavior of CR diffusion on relatively small 
length scales, we use not a concrete GMF model but 
approximate its regular component by a uniform field. 
We measure the momentum distribution /(p) of CRs 
crossing a set of spheres with radii between 1 pc and 
1 kpc centered on the CR source. We verify that in 
the quasi-Gaussian regime the distribution is compati¬ 
ble with a dipole, and compute the anisotropy A of the 
flux of CRs crossing each sphercQ- The statistical error 
of the anisotropy estimate is a factor 1/A^ larger than 
the error of the corresponding flux estimate. Thus even 
for relatively large flux anisotropies, A ~ 0.01, the num¬ 
ber of trajectories and hence the CPU time required is 
increased by a factor 10^. The average total computing 
time used for a single run was about 10 CPU years. To 
further reduce the statistical uncertainty, we have aver¬ 
aged the anisotropy over regions selected according to 
the symmetry of the regular magnetic field Bq: spherical 
symmetry for Bq = 0 and axial symmetry for Bq = const. 

In the left panel of Fig. [TJ we show the flux anisotropy 
A of single sources at various distances as function of 
time for a purely turbulent magnetic field with strength 
Brms = 0.3 pG. We choose the times and length scales 
sufficiently large such that the CR propagation proceeds 
in the quasi-Gaussian regime, t :s> 2Dj and I 3> ^coh- 
For a purely turbulent field, CRs should perform a ran¬ 
dom walk. Thus their density is Gaussian , and as a re¬ 
sult, Eq. ([T]) simplifies (|Shen fc Maolll97ll) to 

A=^.3.3xlO-f^V^V\ (2) 

cT \200-pcJ V2Myr7 ^ ' 

The flux anisotropy A shown in Fig. [T] for sources at 
varying distances is rescaled to the value of the nearest 
source at 50 pc assuming the validity of Eq. Within 
the numerical precision of our results, the obtained values 
for A follow nicely the analytical result. 

Next we add a regular field to the turbulent field. In 
this case, CRs diffuse faster parallel than perpendicu¬ 
lar to the regular field lines and one may wonder, if the 
simple formula (l2|) for the dipole anisotropy still holds. 

^ Note that the anisotropy A = (Fmax — -^min)/(-^max + 
of the CR flux F is connected to the anisotropy (5 = (/max — 
/min)/(/max + /min) of the CR intensity I by 5 = 3A/2. 


Restricting ourselves again to the case t :$> 2n/c^ and 
I ^ ^coh, we confirm numerically that CR propagate 
also with a non-zero regular field in the quasi-Gaussian 
regime. Now the widths ADA of the three-dimensional 
Gaussian describing the CR density n{E,x) are given 
by the corresponding components of the diffusion ten¬ 
sor. Consequently, we expect that the regular field leads 
to strong deviations from a radially symmetric number 
density of CRs but does not influence the anisotropy. In 
the right panel of Fig. [U we show our results for the flux 
anisotropy A in the presence of both a turbulent and 
regular field, choosing as an example the field-strengths 
Bq = 1 /rC and i?rms = 0.3 /rC. Comparing the results 
in the left (without) and the right (with regular field) 
panel, it is clear that the regular field has no impact on 
the anisotropy. 

3. INTERPRETATION OF THE OBSERVED CR 
ANISOTROPY 

In Fig. [21 we s how as band the range of experimental 
data collected bv iDi Sciascio fc: lupoal (|2014ll on the mag¬ 
nitude of the CR dipole anisot ropyci together with data 
from IceCube (blue errorbars) ([Abbasi et al.ll20ll l. The 
anisotropy grows as function of energy until E ^ 2 TeV, 
remains approximately constant in the range 2-20 TeV, 
before it decreases again. Starting from 100 TeV the 
anisotropy increases, with a rate which is consistent with 
one determ ined in the escape model of IGiacinti et al.l 

In the standard diffusion picture, the anisotropy is 
connected to the gradient of the global CR sea density 
and the diffusion coefficient. Both quantities increase 
with energy and, consequently, the energy dependence 
shown in Fig. [2] is difficult to explain. Our results from 
the previous section suggest to connect the plateau in 
the dipole anisotropy between 2-20 TeV to the presence 
of a single nearby source. We determine the relative 
contribution to the total CR intensity from the nearby 
source, fs = I„(E)/ Ent, and from the glob al CR sea 
using the fluxes from IKachelriefi et al.l (|2ni5h . The con¬ 
tribution of the local source to the dipole anisotropy fol¬ 
lows then as 5s = 3fiR/{2cT). This anisotropy is shown 
with red lines in Fig. |2| Combined with T ~ 2Myr 
as age of the source, we can determine its distance as 
i?~ 200pc X (5J(5 X 10-4). 

Since the spread of the experimental data on the dipole 
is considerable, let us consider for illustration the fol¬ 
lowing two extreme situations: First, the true dipole 
may be close to the lower limit of the measured range, 
^ ^ 3 X 10-4. In this case, the plateau extends up 
to ~ 10^4 eV, the maximal energy of CRs accelerated 
in the local source is unrestricted and the source dis¬ 
tance could be as low as 100 pc. The following rise of 
the dipole is then explained by the contribution of the 
gl obal CR sea which we cal culate in the escape model 
of iGiacinti et al.l (|20l4 l2M5il . In the other extreme, the 
true dipole may be close to the upper limit of the mea¬ 
surements, 5 ^ IQ-^. In this case, we have to introduce a 
cutoff exp(—U/Umax) in the energy of CRs accelerated in 
the local source such that their contribution to the total 
CR intensity and thus to the observed dipole decreases 
towards the FAS-Top value. 

We can constrain these extreme choices consider¬ 
ing also the total CR proton flux. In Fig. [^ we 
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Fig. 1.— The rescaled flux anisotropy A of single sources at various distances as function of time. The distance is indicated by the color 
code in kpc. Left panel: turbulent magnetic held Srms = 0.3 pG; right panel: Bq = 1 pG and Brms = 0.3 pG. 



Fig. 2. — Lower and upper limit (green band) on the dipole 
anisotropy and data from IceCube (blue errorbars) compared to 
the contribution from the local source (red, for two values of i?max) 
and from the average CR sea (magenta) as function of energy. 



Fig. 3. — Gontribution from average Galactic (blue line), from the 
local source (red lines, for cutoff i?max = 10^'^ and 10^® eV), and 
old local sources (green line) to the GR proton intensity together 
with experimental data (errorbars). Additionally, the flux from the 
local source is shown by dots for three source ages. 


compare e xperimental resul ts for the pro ton flux from 
CRE AM (lYoon et al.l 120111) . PAMELA (lAdriani et al ' 


20n), KASCADE and KASCADE-Grande (lAoel et al 


20111) . to the proton flux from the local source shown by 


dots and the average global CR sea at low energies (green 
line). Larger maximal values of <5 require a lower value of 
Amax, making the source spectrum more “bumpy”. Such 
a spectral distortion in the total proton flux is avoided 
for values of Amax ^ 10^^ eV, what also ensures that the 
antiproton data from AMS-02 can be explained. 

Let us now comment on the decrease of the measured 
dipole anisotropy at low-energies, < 2 TeV. A natural 
explanation for this decrease is a small off-set of the 
Earth with respect to magnetic field line going through 
the source. The perpendicular extension of the elongated 
volume filled with CRs by the source decreases with en¬ 
ergy. Thus the lower end of the plateau is determined 
by the perpendicular distance of the Earth to the mag¬ 
netic field line going through the source and is therefore 
a free parameter. In Fig. [3l we show the proton flux 
from the local source as a red line using Amin = 10^^ eV 
and Amax = 10^“^ eV and 10^® eV as low- and high-energy 
cutoffs, respectively. 

The difference at low energies between the observed 
CR flux and the one from the local source should be 
contributed by old, local sources. This contribution is 
shown for illustration in Fig. 3 as a green line, obtained 
by subtracting the two fluxes. Fry et al. (2015) showed 
that their age should be T > 14 Myr. Therefore the 
high-energy CRs emitted by these sources escaped al¬ 
ready, what might explain the suppression of their com¬ 
bined CR fluxes above 10^^ eV. Their individual dipoles 
oc fi/T are smaller, should partially cancel in their 
sum, and we therefore neglect their contribution to the 
total dipole amplitude in Fig. [21 Since the number N 
of sources contributing is large, the approximation of a 
smooth source distrib ution and thu s also t he standard 
appro a ch used e.g. bv IStrong et all (I2007D : lEvoli et ^ 
(I2008D : IPutze et al.l (|2011l) mav be justified in this en¬ 
ergy range. The transition between the single and many 
source regimes around A ^ 2 x 10^^ eV may be connected 
to the spectral breaks observed by PAMELA. 

The phase of the first harmonics is almost constant 
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up to 100 TeV energy (|Di Sciascio fc Iud^I2014D . and 
changes thereafter fast. Such a behavior is natural in 
our model: At all energies E < lOOTeV, the anisotropy 
is dominated by local source(s) which are located pref¬ 
erentially along the local magnetic field lines. As result, 
the dipole phase expected does not change its value up 
to lOO TeV, where the Galactic sea of CRs starts to dom¬ 
inate and thus the dipole direction is given by Eq. O- 

The presence of a nearby source dominating the CR 
flux and dipole in the Te V range is natural to ex¬ 
pect (jKachelriefi et al.]l2015f) : The average rate of SN ex¬ 
plosions in the Milky Way is one SN per (0.3-3) x 10"^ yr 
per kpc^. Cosmic rays in the TeV energy range fill a 
100 pc wide, kpc long volume directed along the regular 
GMF lines for a time of a Myr. It is thus likely that one 
SN has exploded during the last Myr within this volume. 

Such a source is consistent with the explanation for the 
deposition of ®°Fe isotopes in the deep ocean crust by the 
passage of an expanding shell of a 2 Myr old supernova 
remnant through the Solar Syste m (|Knie et al.l 119991 
IBenitez et al.ll2Ci02l: iFrv et al.ll2015ll . The distance deter- 
mined bv iFrv et al.l l 201511 e.g. for the case of a 25Mq 
core-collapse supernova is ~ 130 pc. Relaxing some of 
the assumptions like an isotropic emission, this estimate 
may be compatible with our value R ^ 200 pc in case of 
a low dipole amplitude. Moreover, a single local source 
dominating the local CR proton spectrum explains the 
known differences in the slopes between CR protons and 
nuclei and the anomalous energy dependence of the spec - 
tra of positrons and antiprotons (iKachelriefi et al.]l2015ll . 


4. DISCUSSION AND CONCLUSIONS 

The standard approach to CR propag ation using ei¬ 
ther ana lytical (IBerez inskii et al.l 119 9011 or numerical 
methods (fStrong fc Moskalenkolll998HStrong et al.l[2007t 


lEvoli et al.ll200^ assumes a C R source distrib u tion th at 

is smooth in time and spa ce. Blasi fc Amatol (|2012al lbll 
based on earlier work by iLeeT i 1^911 studi^ the fluc¬ 
tuations in the CR density and the dipole induced by 
the stochastic nature of CR sources. Their results for 
the fluctuations diverge in the limit of small source dis¬ 
tances and ages, i?min,7min —>■ 0. These divergences in¬ 
dicate that a statistical description of quantities like the 
CR dipole anisotropy in terms of its ensemble average 
and variance is not adequate, since the results depend 
strongly on the actual properties of the nearest source(s) 


in a given realization. 

The strongly anisotropic diffusion of CRs that is typ¬ 
ical for the GMF parameters favored by the escape 
model enhances these fluctuations: The majority of re¬ 
cent nearby sources is not connected to us by the regular 
field, and their contribution to the local CR intensity is 
therefore suppressed. In contrast, the flux of the single 
(or the few) sources active in the last few million years 
that are located in the volume aligned with the GMF 
lines passing near the Solar system is strongly enhanced. 
These effects prevent the mixing of CRs with energies 1- 
100 TeV of various sources into a “global, average CR 
sea,” and undermine thereby the basic assumption of 
a smooth, global CR distribution inherent in most ap¬ 
proaches to Galactic CR physics. 

We have shown that the dipole anisotropy in the CR 
flux emitted by a single source is independent of the tur¬ 
bulent and regular magnetic field and of the CR energy 
in the quasi-Gaussian regime. In particular, we have 
shown that the simple formula A = R/ {cT) for the dipole 
anisotropy holds also including a regular field. As an ap¬ 
plication we have considered the experimental data on 
the dipole anisotropy. We have argued that the approx¬ 
imately energy-independent plateau in the anisotropy 
around 2-20 TeV can be explained by the presence of 
a nearby CR source. The age and the distance of this 
source are compatible with the one required to explain 
the “anomalies” in the CR intensities of protons, antipro¬ 
tons and positrons. If the source is also responsible for 
the observed ®*^Fe overabundance in the million year old 
ocean crust, its distance should 100-200 pc, favoring thus 
values of the dipole anisotropy close to the lower end of 
the measured values. 

Finally we note that our results can be used also to 
study the dipole anisotropy of the electron flux, con¬ 
straining e.g. the contribution from recent pulsars, if en¬ 
ergy losses are taken into account. 
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